Non-equilibrium Born-Oppenheimer potential energy surfaces for molecular wires 
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We present a method for computing non-equilibrium, current-dependent Born-Oppenheimer po- 
tential energy surfaces for molecular wires. Calculations are performed for polyacetylene wire de- 
scribed by tight-binding model with electron-phonon interactions. We find that dimerization of 
the polyacetylene wire is amplified by electric current fiow. We show that the boundary between 
transparent and opaque black states of the wire is blue-shifted by current. 
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One of the major goals in modern nanotechnology is 
the construction of an electronic circuit in which molec- 
ular systems act as conducting element jlj. Two critical 
issues should be thoroughly investigated to proceed with 
the design of real molecular devices: heating in a molecu- 
lar junction and current induced conformational changes 
of a junction caused by electromigration. By the inelas- 
tic scattering of electrons the kinetic energy of current is 
transfered to the nuclear vibrations and then it is released 
as heat. The second process, electromigration, is the phe- 
nomenon of atom motion due to current induced forces 
and it is related to the transfer of electronic momentum 
to nuclei due to both elastic and inelastic scattering. The 
electromigration can lead to device breakdown at criti- 
cal values of current and therefore current-induced forces 
may present a limitation for the development of molec- 
ular electronic devices. Very recently, controlled electro- 
migration has been used constructively to make a single- 
molecular junction 0. Both issues, heating and electro- 
migration, are ultimately related to understanding how 
Born-Oppenheimer (BO) potential energy surface of a 
molecular wire is responded to electronic current flow. 
If current-dependent BO potential energy surface of a 
molecular device is available, we can compute current- 
induced forces as the corresponding gradients as well as 
we can perform the normal modes analysis to obtain the 
current-dependent vibrations. 

Although the practical importance of a current- 
dependent BO surface can not be underestimated, some 
intricate theoretical questions connected to the difficulty 
in defining a generic variational energy functional for 
non-equilibrium steady state systems |3, 0, ll| can be 
also puzzled out. Theoretical work to date has failed 
to define or even to prove the existence of the BO poten- 
tial energy surface for molecular wires with current flow. 
The principle theoretical difficulty is the treatment of the 
transport problem as an electron transmission (Landauer 
approach ^) which is non- variational from the outset. A 
wire with current is by definition an open quantum sys- 
tem and as long as one works within a non-variational 
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transport theory it is likely not possible to comprise a 
geometrical optimization of the nuclear degrees of free- 
dom with electronic transport calculations within a single 
theoretical approach. There is the opinion that current 
induced forces can not be conservative although va- 
lidity of this claim has been recently questioned j^]- If 
a current induced force is not conservative then there 
does not exist a BO potential energy surface of which 
this force might be derivative. One of the goals of this 
Letter is to show that this is not true and that the Born- 
Oppenheimer potential energy surface can be properly 
defined and efficiently computed in the presence of cur- 
rent flow. 

The Lagrange multiplier based transport theory [To| 
takes its origin in the modern development of steady state 
non-equilibrium statistical mechanics 0, [T3 . . This 
theory is equivalent to a variational formulation where 
the constrained minimization of the expectation value of 
the Hamiltonian on the manifold of the desired current 
J yields the transport properties [§| . The "algorithm" of 
the Lagrange multiplier based transport theory to com- 
pute a current-dependent BO surface can be summarized 
as follows: (1) define the operator of the current J via 
the continuity equation; (2) extend the Hamiltonian H 
by adding the term (— AJ) where A is a Lagrange mul- 
tiplier; (3) compute the BO surface as the expectation 
value (0(A)|iJ|0(A)) where |0(A)) is the "ground state" of 
the extended Hamiltonian H — A J. 

In this Letter, using the Lagrange multiplier based 
transport theory we compute the current-dependent BO 
potential energy surfaces of conjugated polymer, poly- 
acetylene, wire. Polyacetylene is a linear polymer, it 
consists of coupled chain of CH units forming a quasi- 
one-dimcnsional lattice. Each C has four valence elec- 
trons: two of them contribute to the cr-type bonds con- 
necting neighboring carbons along the one-dimensional 
backbone, while the third forms a bond with the the hy- 
drogen. The remaining electron has the symmetry of 
a 2pz orbital and contributes to 7r-bond. In terms of 
the energy band description, the a-bonds form low-lying 
completely filled band, while 7r-bond leads to half-filled 
energy band responsible for electron transport properties. 
The TT-electrons can to a fair approximation be treated 
as a quasi-one dimensional electron gas within the tight- 
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binding approximation. Accordingly, we start our consid- 
eration from the following model electronic Hamiltonian 
for the polymer chain 



2NKx^ . 



(1) 



Here ana creates an electron of spin a on site n. The 
Hamiltonian H is based upon the so called Su-Schrieffer- 
Heeger (SSH) 14| model for the band-structure of trans- 
(CH)„ lattices. The first term in the H gives the energy 
for an electron with spin a to hop between neighboring 
^»z-orbitals. The strength of this term is modulated by 
linear coupling to distortions x in the polymer lattice 
away from evenly spaced lattice positions. Finally, the 
last terms in H gives the harmonic interactions between 
lattice-sites arising from the cr-bonds between neighbor- 
ing lattice atoms. 

We begin the derivation by defining the operator of 
current via continuity equation. The number of electrons 
on the site n is given by the expectation value of the 
operator 



(2) 



By making use of the Heisenberg representation the 
continuity equation can be written as the Heisenberg 
equation-of- motion for the operator iV„: 



Nn - I [Ho, iV„] 



(3) 



Given standard anti-commutation relations between the 
electron creation and annihilation operators, the r.h.s. 
commutator Q can be readily computed. Comparing 
eq.© with the finite difference expression for continuity 
equation Nn = —{jn — jn-i) we obtain the definition of 
the operator of current through site n. By making the 
sum of on-site currents j„ along the wire we define the 
net, macroscopic current through the wire 



(4) 



There are two distinct contributions to the net cur- 
rent. The term to(aJi+io-Q'ncr — alia-O'n+ia-) is the stan- 
dard for molecular wires without electron-phonon cou- 
pling whereas the term {—l)"-ax{al^^j^^ana ~ al^^Un+ia) 
is the phonon-assisted current. 

We assume now that there is a time-independent, con- 
stant current through the wire. A homogeneous current- 
carrying state is the same whether it is introduced by 
reservoirs or by a bulk driving driving field, the current 
enters the problem via a Lagrange multiplier A. To this 
end the Hamiltonian H is modified by adding the term 
which constraints the macroscopic current J: 



Hi 



H -XJ 



(5) 



The Hamiltonian Hj is a quadratic form in fermion cre- 
ation/annihilation operators. This quadratic form can 
be exactly diagonalized by two unitary transformations. 
The Fourier transformation 



ana = — ^ ^exp(ifcn)cfeo-, 

k 

brings the Hamiltonian to the following form 



(6) 



-■k+TTa^ka 



0<fe<27r a 



tk ^k \ I Ck+Tra 1 f~\ 



Cko 



We have introduced here the phonon-unperturbed band 
energy 

efe(A) = 2to(cos(fc) -f Asin(fc)), (8) 

and the phonon-induced gap 

Afe(A, x) = -4ax(sin(fc) - Acos(fc)). (9) 

Both quantities, efe(A) and Afe(A,x) depend upon cur- 
rent via the Lagrange multiplier A and if we let A tend to 
zero we recover the usual zero-current results 14]. Mix- 
ing Cka and Cfc+TTCT operators by the canonical Bogoliubov 
transformation 



Oik+-Ka 
Ctka 



Uk Vk 
Vk -Uk 



^k+ira 
Cka 



(10) 



we obtain new quasipaticles aka and ak+Tva- The co- 
efficients of the Bogoliubov transformation is chosen in 
such a way that the extended Hamiltonian (j?)) should be 
diagonal. It leads to the following expressions for the 
coefficients 



Uk = A/6'(efc)cos(0fc) - \Jl- 9{ek)sm{(j)k) , (11) 



Vk = y/e(7kjsm{(t>k) + y/l~9{ek)cos{(t>k) , (12) 

where ^(e^) is the Heaveside step function and the mixing 
angle (j)k is given by the formula (/)fe = 1/2 arctan(efe/Afe). 
This canonical Boguluibov transformation produces the 
non-interacting Hamiltonian for the quasiparticles aka 
and ak+TTa- 



Hj = 



E E 

Q<k<TT a 

2NKx^ , 



Ek{4 



k+ira'^k+Tra 



"ka 



aka) 



(13) 



with the wave-vector k in the reduced Brillouin zone 
of the system. The band energy, Ek = ^^e^ + A^, 
and the mixing angle (pk in the canonical Bogoliubov 
transformation are functions of the phonon-unperturbed 
band energy, efe(A), and the phonon-induced energy gap, 
Afc(A,x). 

It can be straitforwardly shown that this diagonaliza- 
tion of the Hamiltonian is equivalent to the variational 
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FIG. 1: Born-Oppenheimer potential energy surfaces for dif- 
ferent values of the current: 1=0.00 mA (lower panel), 1=0.25 
mA (middle panel) and 1=0.50 mA (upper panel). The dimer- 
ization coordinate, Xo, corresponds to the minimum on the 
BO surface and it is shown for each values of the current. 



problem 5{H — A J) = with the coefficients of the Bo- 
goliubov transformation Uk and Vk as variational param- 
eters. Current carrying steady states are by definition 
the ground states of extended, current dependent Hamil- 
tonian Hj Jjj . The ground state of the Hamiltonian 
can be defined as the vacuum state for quasiparticle 
ak+TTu and . All physical properties of the system de- 
pend upon the current and can be obtained as the matrix 
element of the corresponding operators over the current 
carrying "ground state" |0(A)). Now the non-equilibrium 
Born-Oppenheimer potential energy surface is defined as 
the expectation value of the original Hamiltonian H 
over the current carrying wave function: 



Ebo{\x) = (0(A)|i/|0(A)) 
= -2 ^ Ek{\x) +2NKx'^ 



(14) 



0<fc<7r 



The band energy Ei^^ the mixing angles 1?^^, and the 
current carrying wave function |0(A)) are not yet in a 
form to be computed as the Lagrange multiplier A is not 
known yet. The additional equation for the Lagrangian 
multiplier A is obtained if the density of an expectation 
value of the net current operator Q over the current 
carrying state of the Hamiltonian is required to yield 
a desired current density /. This sets up the continuity 
equation: 



-(0(A)|J|0(A)) = - (2f?(£fe)-l)(2iosin(fc)cos(2,^fc)+4axcos(fc)sin(20fc))=/ (15) 

0<fc<Tr 

I 



For any given lattice distortion x the nonlinear equation 
(|15|l should be resolved for A. Then the band energy and 
the mixing angle (pk of the canonical Bogoliubov trans- 
formations arc computed as functions of current. 

We now discuss the main results obtained using the 
above model. The numerical calculations are performed 
for the half-occupied conductance band with the box 
length N — 1000. In all our calculations we use the fol- 
lowing model parameters of the SSH model ^^: to = 2.5 
eV, a = 4.1 eV, and K — 21 eV/ A^. This parameteriza- 
tion of the ASH Hamiltonian is typically used to describe 
coherent electron transport in polyacetylene wire p^ . 

In the first calculations, shown in Fig.l, we studied 
the changes of the BO surface as current flows through 
the wire. It is immediately evident from Fig.l that the 
current substantially alters the Born-Oppenheimer po- 
tential energy surface of the polyacetylene wire. The 
minima on the BO surface correspond to the optimal 
non-equilibrium steady state geometries of the wire. The 
dimerization coordinate increases as current flows, i.e. 
the C = C double bonds are shortened and the C — C 



single bonds are elongated by the current induced forces. 
By increasing the current we also increase the barrier 
between the two- fold degenerate ground states, there- 
fore the energy difference between evenly spaced poly- 
acetylene and dimerized polyacetylene is increased as we 
increase the current. Since SSH Hamiltonian is very 
general and applicable in its present formulation to any 
long semiconductor wires we do expect that effect of 
the current-amplified dimerization occurs in other quasi- 
one-dimensional wires as well. The results of our cal- 
culations are qualitatively consistent with the finding of 
Todorov, Hoekstra and Sutton |16|| that current-induced 
forces generally form an alternating patterns in atomic 
wires. The critical value of the current at which the cur- 
rent induced forces break the wire into pieces can be also 
determined. If we assume that the current-induced bond 
distortion larger than 0.2A will be sufficient to rupture 
the wire, the critical value of the current which the poly- 
acetylene wire can sustain is 0.8 mA. 

Finally, we turn our attention towards what physi- 
cal quantities which are accessible experimentally can 
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FIG. 2: The current-dependent valence and conductance 
bands of polyacetylene wire, Ek = ±-^e| + A^. 
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FIG. 3: The transparent-opaque boundary for polyacetylene 
wire as a function of current 



be used to monitor the current-amplified polyacetylene 



dimerization. To this end we compute the response 
of the band structure of polyacetylene to current flow. 
The current-dependent band structure of the first Bril- 
louin zone is plotted in Fig. 2. Polyacetylene's band is 
half occupied therefore the states with negative Ek are 
completely occupied and the states with positive Ek are 
empty. The current flow through the wire does not only 
produce the overall shift of the band energy within the re- 
duced Brillouin zone but also changes the the total band- 
width from 10 eV at zero current to 12.5 eV at J = 0.5 
mA. The band gap between the valence and conductance 
bands is also modulated by current. 

From the band structure calculations the band gap can 
be extracted as a function of the current. Then the band 
gap can be turned into the longest possible wavelength 
of absorbed radiation Xmax and this is plotted in Fig. 3 
for different values of the net current density. Being a 
semiconductor wire, polyacetylene is transparent to the 
photons with wave length longer than Xmax (i-e. with 
the energies within the band gap) but it is opaque to 
photons with wavelength shorter than Xmax (within the 
bandwidth). It is clear from that the boundary between 
transparent and opaque black parts of the spectrum is 
blue-shifted by current flow in the polyacetylene wire. In 
other words, the color of the polyacetylene wire can be 
modified by electric current. 

In this Letter, we have laid a theoretical ground-work 
for computing a Born-Oppenheimer potential energy sur- 
face for molecular wires with direct current. Within the 
Lagrange multiplier based transport theory we presented 
a computational tractable theoretical scheme to compute 
BO of molecular wire junctions. Based upon applications 
of the techniques to tight-binding model with elcctron- 
phonon interactions we predict that (l)the dimerization 
of the polyacetylene wire is amplified by current; (2)the 
transparent-opaque boundary is blue-shifted by current 
flow; (3)the generic mechanism of the semiconductor wire 
breakage is current increased dimerization. 

The author wishes to thank A.Nitzan for valuable dis- 
cussion. 
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